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Dispersed clay particles in mine tailings and soft sediments remain suspended for 
decades, hindering consolidation and challenging effective management of these aqueous 
slurries. Current geotechnical engineering models of self-weight consolidation of tailings 
do not consider microbial contribution to sediment behavior, however, here we show 
that microorganisms indigenous to oil sands tailings change the porewater chemistry and 
accelerate consolidation of oil sands tailings. A companion paper describes the role of 
microbes in alteration of clay chemistry in tailings. Microbial metabolism in mature fine 
tailings (MFT) amended with an organic substrate (hydrolyzed canola meal) produced 
methane (CH 4 ) and carbon dioxide (CO2). Dissolution of biogenic CO2 lowered the pH 
of amended MFT to pH 6.4 vs. unamended MFT (pH 7.7). About 12% more porewater 
was recovered from amended than unamended MFT during 2 months of active microbial 
metabolism, concomitant with consolidation of tailings. The lower pH in amended MFT 
dissolved carbonate minerals, thereby releasing divalent cations including calcium (Ca 2+ ) 
and magnesium (Mg 2+ ) and increasing bicarbonate (HCO3) in porewater. The higher 
concentrations increased the ionic strength of the porewater, in turn reducing the 
thickness of the diffuse double layer (DDL) of clay particles by reducing the surface charge 
potential (repulsive forces) of the clay particles. The combination of these processes 
accelerated consolidation of oil sands tailings. In addition, ebullition of biogenic gases 
created transient physical channels for release of porewater. In contrast, saturating 
the MFT with non-biogenic CO2 had little effect on consolidation. These results have 
significant implications for management and reclamation of oil sands tailings ponds and 
broad importance in anaerobic environments such as contaminated harbors and estuaries 
containing soft sediments rich in clays and organics. 
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INTRODUCTION 

The oil sands in Alberta, Canada are the world's third largest 
proven reserve comprising ~170 billion barrels of recover- 
able bitumen, with production expected to surpass 3.5 million 
barrels day -1 by 2020. Bitumen is extracted from surface- 
mined ores using hot water and hydrocarbon solvent (Schramm 
et al, 2000), generating ~1 million m 3 fluid fine tailings 
day -1 that are deposited and retained in tailings ponds. The 
current total volume (>920 million m 3 ) and surface area 
(~182 km 2 ) of the ponds will continue to increase with resource 
exploitation (Alberta Environment and Sustainable Resource 
Development, 2013) (http://www.aer.ca/rules-and-regulations/ 
directives/tailings-plans-2012). Two major challenges plague the 
oil sands surface mining industry: ( 1 ) recovering water from tail- 
ings for re-use (7.5-10 barrels of water for every barrel of bitumen 
produced), so as to reduce the need for fresh water, and (2) con- 
solidating the tailings solids to decrease the volume of stored 
tailings for subsequent reclamation. 



The tailings deposited into ponds (initially at ~10wt% 
solids) settle by gravity in 2-4 years to form mature fine tail- 
ings (MFT, ~30wt% solids), a colloidal suspension comprising 
slightly alkaline water, silt, clay, unrecovered bitumen (~5wt%) 
and residual solvent (<lwt%) such as naphtha, a mixture of 
C3-C14 aliphatic and monoaromatic hydrocarbons. MFT then 
requires decades to consolidate to >60 wt% solids (Fedorak et al., 
2003; feeravipoolvarn et al., 2009) through close packing of clay 
particles accompanied by interstitial water (porewater) expression 
to form a "cap water" layer. Tailings ponds harbor indigenous 
microbial communities (Penner and Foght, 2010; Siddique et al., 
2011, 2012) that anaerobically biodegrade organic compound, 
such as solvent hydrocarbons, to methane (CH4) plus carbon 
dioxide (C0 2 ) (Siddique et al, 2006, 2007, 2011). Biogenic gas 
production has been associated with accelerated consolidation 
of MFT and recovery of porewater both in situ and in the lab- 
oratory (Fedorak et al, 2003; Bressler et al, 2010). To date, 
only simple physical mechanisms such as creation of transient 
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channels for water transport by biogenic gas bubbles (Brown 
et al., 2013; Voordouw, 2013) have been proposed to explain 
this phenomenon, but this mechanism inadequately describes 
consolidation when gas is formed but ebullition is minimal. 

Because biodensification is slow in situ when supported solely 
by endogenous substrates, we amended MFT with low concentra- 
tions of labile organic substrates to enhance anaerobic microbial 
activity and accelerate the processes ex situ. Here we describe a 
laboratory study conducted using 50-L columns (Figure 1). The 
results reveal biogeochemical reactions that impact the chem- 
istry of three tailings components — porewater, expressed pore- 
water (cap water) and solids — and accelerate consolidation of 
MFT. We describe how microbial metabolism alters the chem- 
istry of porewater that in turn influences consolidation of clay 
particle suspensions. A companion paper (Siddique et al., 2014) 
describes microbe:mineral interactions in the solids fraction. The 
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FIGURE 1 | Experimental 50-L columns used for measuring 
consolidation, porewater recovery and gas release from oil sands 
mature fine tailings (MFT) either amended (A) with an organic carbon 
source or unamended (U). The photograph was taken after 105 d 
incubation at — 20°C. Three sampling ports for collecting MFT or cap water 
(A1-A3; U1-U3) and three instrumentation ports (P1-P3) fitted with 
pressure transducers and a pressure reading unit (PR) for recording in situ 
pore pressure were installed in each column; see Methods for details. The 
columns were sealed under a N2 gas atmosphere. Stainless steel tubing 
connected the headspace to a micro-gas chromatograph (GC) to determine 
the composition of released biogenic gas, and to gas traps (GT) for 
measuring cumulative emitted gas volumes. 



combined results and biogeochemical models are relevant for 
proposed reclamation strategies for oil sands tailings in uplands 
and end pit lakes. 

In addition, the fundamental mechanisms reported here have 
broad implications for developing geotechnical models to predict 
behavior of organic-rich saturated soft sediments (i.e., contami- 
nated harbors, and marine, lacustrine and estuarine sediments), 
predicting contaminant transport from sediments to overlying 
water, and assessing remedial strategies for containing contam- 
inants through in situ sediment capping (Himmelheber et al, 
2011). 

MATERIALS AND METHODS 
PREPARATION OF 50-L COLUMNS 

Mature fine tailings (MFT; Table 1) and cap water were col- 
lected from Mildred Lake Settling Basin (MLSB; 57°4'27"N, 
111°38'19"W UTM) at Syncrude Canada Ltd. in Fort McMurray, 
Alberta, Canada in April 2011, and stored at room tempera- 
ture (~20°C, ca. tailings pond in situ temperature) for a few 
days before using for the experiment. The solids and bitumen 
contents of the MFT were determined using the Dean and 
Stark method (Dean and Stark, 1920) to be ~35 and ~2.5%, 
respectively. The solids content was adjusted to 25wt% solids 
using tailings pond "cap water" immediately before beginning 



Table 1 | Chemical characteristics of mature fine tailings used in 50-L 
column experiment prior to amendment and incubation. 


Analysis (units); number of replicates 


Value 






Bitumen 
Water 
Solids 
Total 


1.5 
74.9 
23.4 
99.8 




Clay <2 |im 
Silt 2-50 |im 
Fine sand 50-250 \im 


25.6 
70.9 
3.5 


BulkpH; n = 2 


7.3 ± 0.03 


E h (mV); n = 3 


-2.0 ± 16 


SOLUBLE CATIONS (mg L" 1 POREWATER) 


Ca 2 +; n = 2 
Mg 2+ ; n = 3 
K+; n = 2 
Na+; n = 2 


18.1 ± 1.1 
8.6 ± 1.2 
24.6 ± 2.7 
737.8 ± 37.6 


SOLUBLE ANIONS (mg L" 1 POREWATER) 


SO 2- ; n = 2 
Ch; n= 3 
HCO3"; n = 3 
PO 2- ; n = 3 


20.6 ± 1.2 
606.0 ± 21.0 
1620 ± 58.5 
BDL C 


Carbonates (wt%) b ; n = 2 


2.4 ±0.1 



Mean values of replicates are shown ± 1 standard deviation. 

a (Dean and Stark, 1920). 

b Calculated on oven dry weight basis. 

C BDL, below detection limit (0.01 mg Lr 1 ). 
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the experiment, to allow this study to be compared to previous 
studies (Young et al., in preparation). For each of two paral- 
lel 50-L columns (described below), 30 L MFT was combined 
with 15.45 L tailings pond cap water by gentle but thorough 
manual mixing under a curtain of N2 gas to achieve a solids 
content of 25 wt% while maintaining anaerobic conditions. For 
the test column, the added cap water contained soluble canola 
meal hydrolysate (preparation and composition described in 
Table 2) sufficient to provide a final concentration of 400 mg 
C L _1 MFT. The cap water used for adjusting MFT in the 
other column was unaltered, serving as a baseline control in 
which only endogenous carbon was available to support microbial 
activity. 

Two cylindrical acrylic columns (30.48 cm outer diameter, 
0.64 cm wall thickness, and 195.6 cm height; Johnston Industrial 
Plastics Ltd., Edmonton, AB, Canada) were sealed at the base 
with acrylic plates and lateral supports (Figure 1). The cylinders 
were equipped with brass sampling ports (U1-U3 for unamended 
and A1-A3 for amended columns) for removal of MFT or water 
samples from the inner portion of the column, and with PX409- 
005G5V (OMEGADYNE Inc.) pressure transducers (P1-P3) and 
an Agilent 34972A Data Acquisition/Data Logger pressure reading 
unit (PR) for measuring pore pressure in situ. 



Table 2 | Elemental composition of hydrolyzed canola meal 3 used as 


the labile organic amendment in the experiment. 




Analysis (units); number of replicates 


Value 




ELEMENTAL ANALYSIS (mg g- l ) b ; n= 1 


Carbon 


28.1 


Nitrogen 


4.3 


Hydrogen 


4.0 


Sulfur 


0.4 




Ca 2+ ; n = 4 


1.6 ± 0.1 


Mg 2+ ; n = A 


3.1 ± 0.1 


K+; n = 4 


15.7 ± 0.2 


Na+; n = 4 


291.3 ± 5.6 




SO 2 ."; n = 2 


2.3 ± 0.1 


CI - ;/) = 2 


366.9 ± 11.2 


HCO§-; n = 2 


213.3 ± 8.9 


P0 2 -;n = 2 


0.9 ± 0.04 



Mean values of replicates are shown ± 1 standard deviation. 
"Canola hydrolysate was prepared from dry canola meal (post-pressing to 
remove oil; provided by Sanimax, Canada) by homogenizing the meal in a mortar 
and adding 10 g of powder to 90 ml of 4% sodium hydroxide. The slurry was pre- 
heated to 60° C for 5min and then stirred at 180 rpm for 24 h at 6CPC. Following 
hydrolysis, the residue was collected by centrifugation at 9000 g and the super- 
natant was collected and adjusted to pH 9 with drops of 1 M HCI. The resulting 
liquid was freeze-dried and stored as a powder at room temperature. 
b Calculated on oven dry weight basis. 

c To assess the contribution of the canola meal amendment to soluble cations 
and anions in MFT porewater and cap water, the hydrolysate was dissolved 
in nanopure water at the same ratio as that used in the amended MFT 
(400mgCL- 1 ). 



The experiment was established on May 6, 201 1 by filling each 
column with 45 L prepared MFT under a plastic curtain flushed 
with a continuous stream of nitrogen (N2) to maintain anaer- 
obic conditions and minimize exposure of MFT to air. In both 
columns, the adjusted MFT formed a single layer of suspended 
clays with no free water at the surface. After filling, the columns 
were sealed under an N2 headspace using acrylic plates fitted 
with gas-tight stainless steel gas ports and tubing (0.15 cm diam.) 
connecting the headspace to a gas trap and a micro-gas chro- 
matograph to determine the volume and composition of emitted 
gas, respectively. The columns were incubated in dim natural 
light at ~20°C. Real time measurements of biogenic gas pro- 
duction and composition, in situ pH, porewater recovery and 
consolidation of tailings were recorded during 213 d incuba- 
tion. Consolidation and porewater recovery were measured as 
the height of the cap water-solid interface (the "mud line") and 
corrected for trapped gas volume as described below. 

At the end of the experiment, cap water and MFT samples 
were collected from the center of the column via sampling ports, 
using 4 mm stainless steel tubing and syringes. Samples were ana- 
lyzed for carbonate minerals, soluble and exchangeable cations, 
and soluble anions, as described below. These samples were also 
subjected to solid phase chemical and molecular microbiological 
analyses described and discussed in the companion manuscript 
(Siddique et al, 2014). 

We were unable to include sterilized MFT as a control con- 
dition for several reasons. First, we have observed previously 
(unpublished results) that it is very difficult to ensure steriliza- 
tion of the thick MFT suspension, even with lengthy, repeated 
autoclaving over three or four successive days. Second, we did not 
wish to alter the structure or chemistry of the tailings by subject- 
ing them to heat, as this could affect consolidation. Finally, we 
did not have access to a gamma source of sufficient capacity for 
irradiative sterilization of ~90 L of MFT. 

SMALL-VOLUME COMPLEMENTARY EXPERIMENTS 

To observe the effect of non-biogenic (exogenous) gas amend- 
ment on consolidation of MFT, four parallel 2-L columns were 
prepared. Each column was filled with 1.7 L MFT (adjusted to 
25% solids) under anaerobic conditions. One column received 
MFT that had been sparged with N2 gas (5.0 Ultra high purity, 
Praxair) at a flow rate of 62.5 mL min -1 for 1 h prior to dispensing 
into the column. MFT for two replicate columns was sparged with 
CO2 (4.0 Anaerobic, Praxair) at a flow rate of 72 mL min -1 for 
1 h until pH stabilized between 6.54 and 6.68, then was placed in 
the columns. MFT for the fourth column did not receive any gas 
sparging (untreated control). All columns were sealed under an 
anaerobic (N2) headspace and incubated undisturbed at ~20°C 
in the dark for 120 d. The columns were monitored regularly for 
expressed porewater recovery (cap water) and solids consolida- 
tion by measuring the height of the mud line and water surface 
and converting to volume. 

To test the effect of different cations, ionic strength and pH on 
the settling behavior of clays particles separately from microbial 
effects, a small-volume, short-term experiment was conducted 
in triplicate in 15-mL conical Falcon tubes (Fisher Scientific, 
Canada). Seven milliliters of MFT was placed in each Falcon 
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tube and amended with 7 mL of one of the following solutions: 
IN NaCl; IN CaCl 2 ; IN FeCl 2 ; 0.2M ammonium oxalate buffer 
(pH 3); 6M HQ (pH 1); or Nanopure water (Barnstead; Thermo 
Scientific). The tubes were shaken on a reciprocating shaker for 
2 h and then placed upright on the bench to consolidate without 
disturbance. The mud line height in each tube was recorded after 
30 min and again after 24 h. 

ANALYTICAL METHODS AND CALCULATIONS 

We define here two water compartments analyzed in our exper- 
iments: (1) "porewater" is interstitial water associated with the 
colloidal MFT suspension, and is a component of the MFT; (2) 
"cap water" is porewater that has been expressed from the MFT 
and has collected above the mud line during consolidation. The 
term "mud line" used here is defined as the interface between cap 
water and MFT. 

Porewater recovery as cap water, and MFT consolidation 

The height of the mud line was used to measure solids consol- 
idation and calculate porewater recovery. The volumes of cap 
water (V cw ) and MFT (solids plus porewater plus gas below the 
mud line; Vmft) were determined, respectively, by measuring the 
height of the cap water level above the mud line and height of the 
mud line in each column. The measured heights were converted 
to volumes using the following equation: 

V = Itr 2 /! mu d line 

where, V = volume of cap water or MFT, 

r = inner radius of the column, 

b-mud line = height of the cap water level or mud line. 
Then, water recovery (WR) defined as the cap water volume 
(above the mud line) was calculated as a percentage of the initial 
volume of porewater in MFT: 

WR = (V cw /V in ) x 100 

where, WR = water recovery (%), 

V cw = measured volume of cap water, 

Vin = initial volume of porewater MFT. 
Consolidation of MFT is defined as the MFT volume (below the 
mud line) calculated as a percentage of the initial volume of MFT 

Consolidation (%) = [(V in - V M FT)/V in ] x 100 

where, Vmft = measured volume of MFT and Vi n = initial 
volume of MFT. 

Estimation of MFT consolidation in 50-L columns by deducting 
volume of entrapped gas bubbles 

Estimated consolidation of MFT in unamended and amended 
columns was calculated using the following equation: 

Consolidation (%) = [(V in - V M F T _ d )/V in ] x 100, 

where, Vi n is the initial volume of MFT (45 L), and Vmft-cI 
is volume of MFT without gas at different densities that are 



calculated as: 

VMFT-d = (niMFT-In — ni wr )/D 

where, m^FT-in is the initial mass of MFT (54.9 kg); m wr is the 
mass of released water. D is the density of the MFT on different 
days, determined from the mean weight of five replicate 10-mL 
subsamples of MFT. 

Volume and composition of gas emitted by ebullition 

The volume of gas emitted from 50-L columns was measured 
using a NaCl-citric acid brine gas trap (Boone, 1982). The com- 
position of produced gas for CH4 and CO2 was determined 
using a Varian CP-2003 portable Micro-GC (gas chromatograph) 
equipped with a GMT-2HP moisture trap and thermal conduc- 
tivity detector (TCD). The minimum detection limit of this GC 
was 50 ppm. 

pH and redox potential (Ef,) 

The pH was determined in situ at intervals during incubation 
by inserting a pH meter (Hach H170multi) fitted with an ISFET 
pH stainless steel micro probe (PHW17-SS) through side ports 
(Figure 1). For E n measurements, MFT samples were collected 
from ports 2 and 3 under a continuous flow of N 2 and immedi- 
ately subjected to Eh measurement using a TX100 pH/mV-meter 
fitted with redox-sensitive electrode (Sensorex 651CD-ORP). 

Soluble cations and anions 

MFT samples removed from the columns were centrifuged 
(3075 g for 1 h) in a Sorvall RC 5B Superspeed centrifuge 
to collect the porewater. The solids fraction was analyzed as 
described below. For the soluble cations calcium (Ca 2+ ), mag- 
nesium (Mg 2+ ), sodium (Na + ) and potassium (K + ), porewater 
was filtered and diluted with 1% HNO3 (trace metal grade) and 
analyzed using Inductively Coupled Plasma Mass Spectrometry 
(ICP-MS; Perkin Elmer SCIEX ELAN 9000) with appropriate 
internal and external standards. Chloride (Cl~), nitrate (NO^) 
and sulfate (SO 2 ^) soluble anion concentrations were deter- 
mined using an ion chromatograph (Dionex DX 600) equipped 
with a 4 mm analytical column (AS9-HC). Dissolved bicarbon- 
ates (HCO3") were determined using the methyl orange indi- 
cator method (US EPA, 1974) and Smartchem Discrete Wet 
Chemistry Analyzer 200 (Westco Scientific) at 550 nm. The 
PO^ - concentration was determined using the colorimetric 
Molybdenum Blue method (US EPA, 1993) using the SmartChem 
atX880nm. 

Moisture correction factor 

Data from chemical analyses of the solid phase (carbonates and 
exchangeable cations) were corrected for moisture and expressed 
on a solid dry weight basis. Moisture content in the samples was 
determined by drying in an oven at 105°C to constant weight 
(ISO 11465, 1993). The moisture content was obtained in weight 
percent and the corresponding correction factor was developed 
using the following equation: 

Moisture correction factor (MCF) = (100 + moisture %)/100 
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Exchangeable cations 

Five grams of MFT were collected in 50-mL polypropylene cen- 
trifuge tubes (Corning Inc.) and 20 mL of 0.1M BaCl 2 /0.1M 
NH4CI solution was added to displace the exchangeable cations 
(Sparks, 1996; Carter and Gregorich, 2008; Rayment and 
Lyons, 2012). The samples were shaken for 2h and cen- 
trifuged (3075g for 1 h) in a Sorvall RC 5B Superspeed cen- 
trifuge. The supernatant was diluted with 1% HNO3 (trace 
metal grade) and analyzed using ICP-MS. In a few cases where 
the samples were not analyzed immediately after centrifuga- 
tion, the supernatant was preserved by acidification to pH 
1-2 using HNO3. The concentrations of exchangeable cations 
were calculated by subtracting concentrations of soluble cations 
(So et al., 2006). The results were expressed on dry weight 
basis. 

Ionic Strength (I) and Diffuse Double Layer (DDL) 

Ionic strength of the MFT porewater was calculated using the 
following formula (Essington, 2004): 

I = 1/2 £(M,Z?) 

< 

where, M; is the molar concentration of charged species (z) with 
charge Z, in the porewater and the sum includes all charged 
species measured in the porewater. 

The thickness of the DDL was calculated using the following 
formula (Essington, 2004): 

k -1 = (3.042 x 10~ 8 )/Zf°- 5 

where, k _1 is the thickness (cm) of the DDL, Z is the average 
mean charge of the counterions (exchangeable cations), and / is 
the ionic strength of the porewater. Because MFT is a polyionic 
system having ions of different charges, the average mean charge 
of the exchangeable cations was calculated using the following 
equation: 

Z ave = ^ f CfZj/ ^ ^ Cj 

where, Qis the concentration of ions on the exchangeable surface 
and Z, is the charge of the species (i). 

Carbonate minerals 

The total carbonates in the MFT solid phase were determined 
by acid digestion (Pansu and Gautheyrou, 2006). The MFT 
(~10 mL) was centrifuged to separate the solid phase which was 
washed twice with 30 mL methanol to remove any soluble HCO^ . 
The MFT solid phase (~2 g) was placed in a serum bottle capped 
with a butyl rubber stopper and the headspace in the bottle was 
flushed with N2 at atmospheric pressure. The carbonates were 
dissolved by adding 20 mL of 1 M HC1 to the bottles using a 
syringe. The bottle was gently shaken for 2 h and the amount of 
C0 2 released was measured using GC-TCD (HP 5890). After C0 2 
measurement, the contents of the bottle were analyzed for Ca and 
Mg using an Atomic Absorption Spectrometer (AAS; SpectrAA 
880). These concentrations were used to calculate CaC03 and 
CaMg(C0 3 ) 2 in the MFT. 



RESULTS 

We investigated contributions of microbial activity (methano- 
genesis) to consolidation in a series of experiments conducted 
at different scales from 2-L to 50-L column studies. The results 
reveal that methanogenesis is important for rapid consolidation 
of MFT, an effect that is superimposed on tailings self- weight con- 
solidation. Acetate, a well-known substrate for methanogenesis, 
resulted in greater porewater recovery and solids consolidation 
when added to MFT as compared to unamended MFT (Figure 2; 
constructed using data retrieved from Bressler et al, 2010). 
The same phenomenon has also been observed in our previous 
experiments (Arkell et al., in review; Young et al., in prepara- 
tion). The increase in total volume in the acetate-amended MFT 
shown in Figure 2A, left panel is due to biogas entrained in the 
MFT solids. However, specifically inhibiting methanogenesis by 
adding sodium 2-bromoethanesulfonate (BES) resulted in mini- 
mal consolidation, presumably equivalent to self-weight settling, 
minimal porewater recovery and no visible pockets of entrained 
gas (Figure 2). Thus, it appears that active methanogenesis is a 
prerequisite for biodensification of MFT incubated under these 
conditions. 




Porewater 
recovery 



0 20 40 60 

Day 



FIGURE 2 I Effect of adding 10 mM sodium 2-bromoethanesulfonate 
(BES), a specific inhibitor of methanogens, on dewatering of 1.6 L MFT 
in 2-L columns incubated undisturbed for >60 d at ambient 
temperature (~22°C). (A) Photographs of porewater recovery and total 
volume of MFT solids and porewater, including trapped biogenic gases 
(white arrows) in parallel columns containing MFT amended with acetate 
(140 mg L _1 MFT) or with acetate plus BES (few, tiny bubbles; not visible), 
or unamended (no bubbles; not shown). (B) Quantitation of porewater 
recovery from MFT amended with acetate, or with acetate plus BES, or 
unamended, illustrating the magnitude of the microbial effect (Figure is 
constructed using the data retrieved from Bressler et al. 2010). 
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To further demonstrate the contribution of biogenic vs. abi- 
otic gases to porewater recovery and consolidation, small-scale 
(2-L column) experiments were performed in which unamended 
MFT was sparged with CO2 to simulate biogenic CO2 pro- 
duction, or with N2 gas to simulate the effect of insoluble gas 
ebullition (by creating transient physical channels for expression 
of porewater). Only marginal differences in porewater recovery 
and tailings consolidation were observed compared to non- 
sparged MFT (Figure 3). The small proportion and slow rates 
of consolidation and porewater recovery (compared to acetate- 
amended MFT; Figure 2) presumably are due to methanogenic 
activity supported by endogenous organic substrates 
in MFT. 

BIOGENIC GAS PRODUCTION AND DECREASED pH AFFECT 
POREWATER RECOVERY FROM MFT IN 50-L COLUMNS 

In previous small-volume experiments we used complex or 
simple fermentable organic amendments to accelerate dewa- 
tering and consolidation of MFT (Bressler et al., 2010; Arkell 
et al., in review; Young et al., in preparation), and determined 



that a low concentration of hydrolyzed canola meal (400 mg 
C L" 1 MFT) was sufficient to promote tailings biodensifica- 
tion under methanogenic conditions. In the current experiment, 
it is apparent that methanogenic conditions were established 
in the 50-L columns, as anaerobic metabolism of the soluble 
canola hydrolysate resulted in emission of ~2.8 L CH4 and ~0.5 L 
CO2 (CFLirCCh ratio of 5.6:1) to the headspace by 75 d, with 
minor gas production from unamended MFT (Figure 4A). 
Additional biogenic gas was trapped within amended MFT 
(Figure SI Panel B), contributing to the total apparent volume 
of MFT. Stoichiometric calculation of theoretical gas produc- 
tion from canola, based on the Symons and Buswell equation 
(Roberts, 2002), predicted a CH^CC^ ratio of ~3:1, but sub- 
stantially less CO2 than predicted was measured in the headspace 
(Figure 4A). The deficit in emitted CO2 is explained by dis- 
solution of biogenic CO2 in porewater, which decreased the 
pH of amended MFT from 7.7 to 6.4 by 147 d, whereas the 
pH of unamended MFT was essentially unchanged in >200d 
(Figure 4B). We initially observed an increase in pH in the 
amended MFT during the first 35 d, after which it significantly 
decreased. During active gas biogenesis in the first 2 months, 
12% more porewater was released from amended MFT than 
unamended (Figure 4C), most of which coincided with pH 
reduction, a major ebullition event at 50 d, and rapid consol- 
idation of amended MFT. Ebullition of biogenic CH4 created 
transient channels for escape of porewater from MFT, observed 
particularly near the mud line (MFT-cap water interface; Figure 
SI). Other experiments (unpublished) confirmed that transient 
water- and gas-filled channels can form in MFT to allow pas- 
sage of porewater to the surface during this phase. Gas bubbles 
remaining in the amended MFT did not affect subsequent consol- 
idation behavior (Figure S2), but did contribute to the apparent 
total solids volume; complete degassing theoretically would have 
yielded ~15 vol% more consolidation in the amended column 
(Figure 4D). 

pH AFFECTS DISSOLUTION OF CARBONATE MINERALS AND 
POREWATER CHEMISTRY 

The decreased pH in amended MFT due to dissolution of 
biogenic CO2 increased dissolution of carbonate minerals in 
MFT (Figure 5A): at 213 d, the carbonate mineral content in 
amended MFT (1.7 wt%) was significantly lower than una- 
mended MFT (2.7 wt%). Carbonate dissolution increased sol- 
uble cation concentrations in porewater and cap water: Ca 2+ 
and Mg 2+ concentrations were greater in amended MFT (~32 
and 19 mg L , respectively) than unamended MFT (~17 and 
10 mg L _1 , respectively; Figure 5B). In contrast, soluble K + 
decreased in amended MFT (Figure 5B) whereas the solu- 
ble Na + concentration remained unaffected in amended MFT 
porewater (Figure 5C). A very small concentration of solu- 
ble phosphate (PO^ - ) was detected only in amended MFT 
porewater (0.36 mg L _1 ) with low concentrations of sulfate 
(SO 2 " , <5mg IT 1 ) (Figure 5D). Amended MFT had slightly 
higher porewater concentrations of soluble HCO^ (~1500mg 
L~ 1 ) than unamended MFT (~ 1 380 mg L~ 1 ) but Cl _ concentra- 
tions were comparable in both amended and unamended MFTs 
(Figure 5E). 




0 25 50 75 100 125 




0 25 50 75 100 125 

Day 

FIGURE 3 | Effect of abiotic gas purging on MFT consolidation (A) 
and porewater recovery (B). Four replicate 2-L columns containing 1.7 L 
unamended MFT were purged with CO2 (n = 2; error bars, where visible, 
show the standard error) or N2 (n= 1) or no gas purging (n= 1) prior to 
dispensing MFT into columns. Consolidation was calculated from the 
volume of solids below the mud line as a fraction of the total volume of 
MFT plus water plus any trapped gas. Porewater recovery was calculated 
using the recovered porewater volume above the mud line (i.e., cap water 
volume) as a fraction of the initial total volume of porewater in MFT. 
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FIGURE 4 | Physico-chemical changes to unamended (U) and amended 
(A) MFT in 50-L columns. (A) Cumulative emitted methane (CH4) and 
carbon dioxide (CO2) in headspace (not accounting for gases trapped or 
dissolved in MFT); (B) bulk pH of MFT measured at ports 2 and 3; (C) 
Proportion of initial porewater recovered above the mud line as cap water 
(i.e., expressed porewater). The sudden increase in cap water in the 



amended column corresponded to a large ebullition event at ~50 days 
recorded in A; (D) Consolidation of solids, corrected for trapped gas volumes 
by subtracting the estimated volume of trapped gas from the measured total 
MFT volume. Initial and final densities of unamended (D = 1 .22 and 1.30 g 
mL _1 , respectively) and amended MFT (D = 1 .22 and 1.49 g ml_~ 1 , 
respectively) were used to calculate consolidation; see Methods for details. 



EXCHANGEABLE CATIONS. IONIC STRENGTH AND DIFFUSE DOUBLE 
LAYER (DDL) 

Although dissolution of carbonates increased Ca 2+ and Mg 2+ 
concentrations in porewater of amended MFT, no significant 
change in the composition of exchangeable cations on clay sur- 
faces was observed in MFT (Figure 5F). However, increased 
concentrations of soluble cations (Ca 2+ and Mg 2+ ) and anions 
(HCO3") increased the ionic strength (I) of the porewater 
(Figure 6A). The greater I (~0.055mol L _1 ) was calculated 
for the porewater from amended MFT vs. unamended MFT 
(~0.04 mol LT 1 ). The J of a solution has a profound effect 
on the DDL of clay particles. Increased I in the porewater of 
amended MFT decreased the DDL thickness (~10 x 1CP 8 cm) of 
clay particles in amended MFT as compared to the DDL thick- 
ness (~18 x 10~ 8 cm) calculated for clay particles in unamended 
MFT (Figure 6B). 

DISCUSSION 

Dewatering of tailings is a major challenge faced by the 
surface-mining oil sands industry to comply with govern- 
ment directives for managing the ever-growing inventory 



of oil sands tailings (http://www.aer.ca/rules-and-regulations/ 
directives/directive- 074). Any process that increases the settling 
of clay particles in oil sands tailings ponds will enhance consol- 
idation of tailings to reduce the inventory volume and enable 
reclamation, and will increase recovery of porewater from MFT 
for re-use and reduction of fresh water demand; even marginal 
increases in consolidation translate into enormous volumes, given 
the scale of existing tailings ponds (>920 million m 3 ). Thus, 
microbially-mediated consolidation and dewatering (biodensifi- 
cation) can mitigate this problem, whether supported by endoge- 
nous substrates in situ (Fedorak et al., 2003) or by amendment 
with organic substrates in an engineered process ex situ. In this 
paper, we have focused on how microbial activities alter porewa- 
ter chemistry to enhance tailings consolidation and dewatering; 
in the companion paper (Siddique et al, 2014), we consider the 
effect on tailings solids. 

The mechanism of microbially-mediated consolidation (bio- 
densification) of oil sands tailings is multi-faceted, embracing 
different biogeochemical processes occurring simultaneously in 
the MFT (Figure 7). Stimulating the activity of indigenous anaer- 
obes by amending MFT with labile organic carbon (hydrolyzed 
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FIGURE 5 | Concentrations of carbonate minerals and major cations 
and anions in unamended (U) and amended (A) 50-L columns after 
213 d incubation. Labels 1, 2, and 3 refer to samples collected from 
ports on the columns (see Figure 1); Port 1 accessed cap water only by 
213 d; ports 2 and 3 were below the mud line and accessed MFT. Bars 
represent the mean from analyses of duplicate or triplicate samples 



taken from each port and error bars, where visible, represent 1 standard 
deviation. (A) Carbonate mineral content of MFT samples (n = 3). (B,C) 
Concentrations of major soluble cations in cap water and interstitial 
porewater (n = 2). (D,E) Concentrations of major soluble anions in cap 
water and interstitial porewater (n=2). (F) Concentrations of 
exchangeable cations in MFT (n = 2). 



canola meal in this study, acetate in previous studies; Fedorak 
et al., 2003; Arkell et al., in review) enhanced biogenic CH4 and 
CO2 production. Gas production resulted in ebullition of bub- 
bles dominated by CH4 (due to poor solubility of CH4 in water), 
creating transient channels for escape of pressurized porewa- 
ter, particularly in MFT near the mud line. Methanogenesis in 
MFT might be responsible for the observed initial pH increase 
as acetate and CO2 were consumed (Fotidis et al., 2013) but the 
dissolution of entrapped CO2 reduced porewater pH, thereby 
dissolving carbonate minerals and releasing divalent cations 
(Pathway I in Figure 7). Our results are very similar to the 
findings of Wersin et al. (2011) who studied biogeochemical 
processes in a clay formation in situ. They observed carbonate 



dissolution, high pCO>2 and alkalinity (HCO^ and COj~) and 
a decrease in pH from 7.7 to 6.8 during anaerobic degradation 
of an organic substrate under sulfate-reducing and methanogenic 
conditions. Similar results were also reported by Schlegel et al. 
(2011) studying fluid geochemistry and microbiology of multi- 
ple organic-rich reservoirs in the Illinois Basin, USA. Our results 
are also supported by other studies where the influence of benthic 
bacterial activity on carbonate mineral dissolution in marine sed- 
iments has been described (Moulin et al., 1985) and dissolution 
of marine carbonate minerals by acidification due to increased 
atmospheric CO2 has been comprehensively reviewed (Morse 
et al., 2007). Solubility of CaC03 increases 20-fold between pH 
8 and pH 6 (Chou et al., 1989) and siderite can be dissolved 
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FIGURE 6 | (A) Calculated ionic strength (/) of the porewater and (B) 
thickness of diffuse double layer (DDL) of clay particles in unamended (U) 
and amended (A) MFT in the 50-L columns at ports 2 and 3 (Figure 1). Ionic 
strength was calculated using the data presented in Figure 5, the DDL 
thickness was calculated using the ionic strength, and charges of 
exchangeable cations were calculated from data shown in Figure 5F. Bars 
represent the mean values from analyses of duplicate samples taken from 
each port and error bars, where visible, represent 1 standard deviation. 



in anaerobic aqueous environments at pH~7 (Jensen et al., 
2002). 

Dissolution of MFT carbonate minerals (presumably cal- 
cite/dolomite) increased Ca 2+ and Mg 2+ concentrations in 
amended porewater, whereas K + decreased in amended MFT, 
which might be attributed to reducing conditions prevailing in 
amended MFT (Siddique et al., 2014) facilitating K + fixation in 
reduced 2:1 phyllosilicates (Eslinger et al., 1979; Shen and Stucki, 
1994). Phyllosilicate analysis of MFT showed that fine clay (0.08- 
0.2 (xm) comprised ~45% kaolinite-smectite (94% kaolinite) and 
55% illite-smectite (93% illite) group minerals whereas clay (0.2- 
2 |xm) consisted of 62% kaolinite-smectite (100% kaolinite) and 
36% illite-smectite (96% illite) groups (Siddique et al, 2014). 
Among anions, a relatively higher concentration of HCO^ in 
the porewater of amended MFT was observed, presumably due 
to dissolution of biogenic CO2 in porewater and/or dissolution 
of carbonate minerals (Morse et al., 2007). The change in pore- 
water chemistry of amended MFT due to microbial metabolism 
did not alter the cationic composition on clay exchanger sur- 
faces despite significant increases in divalent cation concentra- 
tions (Ca 2+ and Mg 2+ ) in the porewater. However, the divalent 
cation concentrations were dwarfed by the Na + concentration, 
which dominated the porewater composition. Another possible 



reason for the apparent lack of cation exchange could be tech- 
nical, since no method has been optimized yet to determine 
the exchangeable cations in oil sands tailings. Although we can- 
not exclude cation exchange processes based on our results, we 
do not consider cation exchange to be a significant or essen- 
tial process in MFT consolidation (Pathway I, Figure 7, broken 
lines). 

The phenomenon of clay dispersion in aqueous systems is 
understood using the concept of diffuse double layer (DDL), 
the thickness of which is governed by the charge potential on 
the clay surface (one electric layer, generally with net negative 
charge on clay surface) and valence of the counterions in solu- 
tion (second electric layer, positive charge contributed by cations) 
neutralizing the surface charge potential of clay. Decreasing thick- 
ness of DDL increases clay flocculation. In our study, the greater 
ionic strength (I) of porewater from amended MFT, attributed 
to increased concentrations of soluble ions (Ca 2+ , Mg 2+ , and 
HCO3"), would decrease the thickness of the DDL of clay parti- 
cles (Figure 6) and enhance consolidation of tailings. Electrostatic 
repulsion between the clay particles is reduced at higher I due to 
ion shielding of the clay surface charges, decreasing thickness of 
the DDL. Greater consolidation of amended MFT might also be 
attributed to dominance of kaolinite phyllosilicates in the MFT, 
in addition to greater I. The charges on kaolinite are largely pH- 
dependent, and the observed lower pH of MFT would decrease 
the net negative charges on the clay surfaces. Nasser and James 
(2006) examined the settling behavior of kaolinite in aqueous 
media in response to J and found that, at higher I, particles settled 
in flocculated form. Similar results were reported by Mietta et al. 
(2009) for flocculation of a kaolinite suspension as function of 
pH and salt concentration, and they found that clay floe size was 
3-fold greater at pH 4 than at pH 7. Our small-scale, short-term 
(24-h) abiotic experiment conducted with unamended MFT con- 
firmed that adding different salts or decreasing bulk pH by adding 
acid markedly increased consolidation and porewater recovery 
from tailings (Figure S3). Greater consolidation occurred in MFT 
amended with FeCl2 compared to untreated MFT, followed by 
the effects of CaCl2 and then NaCl; maximum consolidation was 
achieved using 6N HC1 (pH 1). These results support our pro- 
posed biogeochemical model (Pathway I, Figure 7), in that: (1) 
increased I of the MFT decreased DDL of clay particles; (2) diva- 
lent cations can exchange monovalent cations from exchanger 
surfaces (although not observed in current study, this was a signif- 
icant effect in an acetate-amended 2-L MFT column study; Arkell 
et al., in review), thereby decreasing DDL; and (3) lower MFT pH 
dissolved carbonate minerals and changed pH-dependent charges 
on clays, leading to consolidation of MFT. 

The results of this study are important for the management 
of oil sands tailings. Microbially-enhanced tailings consolidation 
apparently is occurring in situ (Fedorak et al., 2003), likely sup- 
ported by endogenous substrates such as light hydrocarbons used 
in bitumen extraction, where concomitant recovery of porewater 
can decrease on-site demand for freshwater for oil sand ore pro- 
cessing. Biodensification may also facilitate dry landscape recla- 
mation of MFT by complementing or possibly foregoing current 
chemical and physical management practices like chemical floc- 
culant addition and energy-intensive centrifugation processes. In 
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FIGURE 7 | Proposed model for microbially-mediated geochemical 
pathways of clay consolidation. This study focuses on Pathway I; the 
companion paper (Siddique et al., 2014) focuses on Pathway II. Dissolution of 
biogenic CO2 in tailings decreases porewater pH by producing hydronium ion 
(H30+). Increased H3O+ dissolves carbonate minerals such as CaMg(C03)2 
and releases calcium (Ca 2+ ) and magnesium (Mg 2+ ). Dissolution of CO2 
and/or carbonate minerals increases bicarbonate (HCO^T) ions in porewater. 



Liberated Ca 2+ and Mg 2+ plus HCO3 increase ionic strength (/) of porewater, 
thus reducing the diffuse double layer (DDL) of clay particles and facilitating 
their consolidation. Ca 2+ and Mg 2+ can replace some dominant 
exchangeable Na from clay exchanger surfaces via the cation exchange 
process, reducing DDL of clay particles. Competing reactions are shown by 
solid arrows while dashed arrows indicate pathways not considered 
significant in our study. 



addition to oil sands tailings, it would be prudent to consider 
the roles of indigenous microbes in the consolidation of other 
organic-rich soft sediments, such as those in contaminated har- 
bors and river deltas, and to incorporate biological factors into 
consolidation models. 
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